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文章 编号 :1005-3085(2011)03-0343-11 


一 类 自 催 化 反应 系统 的 Hopf 分 歧 与 稳定 性 * 


马 晓 丽 
(西安 工业 大 学 数理 系 ， 西 安 710032) 
摘 要 : 本 文 考虑 一 类 带 饱 和 项 的 自 催化 反应 系统 . 我 们 首先 讨论 了 党 微分 系统 Hopf 分 歧 的 存在 性 ， 得 
到 了 渐 近 稳定 的 周期 解 ， 其 次 讨论 了 具有 扩散 项 的 偏 微分 系统 ， 在 扩散 系数 满足 一 定 的 条 件 下 ， 
得 到 了 次 临界 的 Hopf 分 歧 的 存在 性 ， 并 且 利用 中 心 流 形 约 化 方法 ， 判 断 出 由 该 Hopf 分 歧 产生 的 
空间 齐 次 的 周期 解 是 渐 近 稳定 的 . 最后， 借助 Matlab 软件 形象 地 验证 和 刻画 了 文中 的 结论 . 
关键 词 : 自 催化 ， Hopf 分 歧 ， 稳 定性， 扩散 项 
分 类 号 : AMS(2000) 35K57 中 图 分 类 号 : 0175.26 文献 标识 码 : A 


1 引言 


本 文 在 Neumman 边界 条 件 下 考虑 一 类 带 饱 和 项 的 自 催化 反应 系统 


us — dy Au = Aww — Ae zEeQ, t>0, 

vt — dgAv = 1 — uv, rE, t>0, o 
1 

Ou = O,v = 0, zE OQ, t>0, 


u(x,0) = uo(z) > 0, #0, v(z,0)=u(z)>0, #0, reEQ, 


EPO R(n > 1) 中 具有 光滑 边界 的 有 界 开 区 域 ，%, v 分 别 代表 两 种 反应 物 的 浓度 ，d1, d2 为 
扩散 系数 ， 和 ,有 为 代表 反应 率 的 系数 ， 并 且 和 , k > 0. 关于 系统 (1) 的 建立 及 相关 的 化 学 背 
景 ， 有 兴趣 的 读者 可 参见 文献 [1,2]. 

对 于 系统 (1) 的 研究 目前 已 有 不 少 工 作 . 在 Dirichlet 边界 条 件 下 ， 文 献 [3] 利用 分 歧 理 论 给 
出 了 正平 衡 态 解 的 存在 性 ， 并 刻画 了 解 的 长 时 行为 ; 文献 中考 虑 了 参数 入 或 上 充分 大 的 情形 ， 
利用 正则 摄 动 和 奇异 摄 动 理论 ， 给 出 了 正平 衡 态 解 的 多 重 性 和 惟一 性 结果 . 在 Neumman 边界 
条 件 下 ， 文 献 55] 利用 不 动 点 指标 理论 给 出 了 非常 数 正 平衡 态 解 的 存在 性 和 不 存在 性 ， 所 得 结 
果 表 明 扩 散 系 数 及 反应 器 的 大 小 在 空间 斑 图 的 形成 中 起 很 大 作用 . 


类 扩散 的 捕食 - 食 饵 模型 的 Hopf 分 歧 ， 给 出 了 空间 非 齐 次 的 周期 解 的 存在 性 及 稳定 性 . 关 
于 Hopf 分 战 的 研究 ， 有 兴趣 的 读者 可 参见 文献 [7,8]. 

大 量 的 文献 表明 ， 对 系统 (1) Hopf 分 歧 的 研究 目前 并 不 多 .鉴于 此 ， 本 文 主要 考虑 自 催化 
反应 系统 (1) Hopf 分 歧 的 存在 性 、 方 向 及 稳定 性 . 同时 ， 借 助 Matlab 软件 ， 从 数值 模拟 的 角 
度 验证 并 加 强 了 所 得 结论 . 


收 稿 日 期 : 2009-09-08. 作者 简介 : 马 晓 丽 (1972 年 5 月 生 )， 女 ， 讲 师 . 研究 方向 ， 数 理 金融 与 微分 方程 . 
* 基 金 项 目 : 国家 自然 科学 基金 (10971124). 
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当 0 < 入 <k-! 时 ， 系 统 (1) 存在 惟一 正常 数 平衡 解 





入 | 


w= (ph 


以 下 均 假 设 0 < 入 < kL. 


2 “” 常 微分 系统 的 Hopf 分 歧 及 稳定 性 


考虑 动力 系统 
u 
T Auv EET t>0, 6) 
dv 
dt => 1 bag UU, t > 0 
经 计算 ， 系 统 (2) 在 (ur, v*) 处 的 雅 可 比 和 矩阵 为 
2 
Te ( XR A 735 | is 
_iAMk a 
入 1 一 入 天 


设 其 特征 方程 为 12 一 Th 十 D =0， 其 中 


A(X2k3 — 2Ak? +k — 1) 
1—Ak ， 





T := trJ = D := detJ = A(1 — Ak) > 0. 


Ad f(A) = A?k —2Ak? +k- 1， 由 求 根 公式 解 得 F(A) = 0 的 两 根 为 
do =k (Vk-1), M = 有 (VE 十 1 


4O<k<1MN, f(A) < 0， 即 有 trJ < 0， 因 此 (w*,v*) 局 部 渐 近 稳定 ， 

定理 1 车 kk > 1， 则 有 以 下 结论 成 立 : 

(i) 当 0 < 入 < 各 时 ，(w*,v*) 不 稳定 ， 当 Xo < 入 <k7! 时 ，(w*,v*) 局 部 渐 近 稳定 ; 

(i) 当 和 = 和 0 时 ， 系 统 (2) 在 平衡 点 (w*,v*) 处 产生 Hopf 分 歧 ， 且 该 Hopf 分 支 为 次 临界 方 
向 ， 并 且 周 期 闭 轨 渐 近 稳 定 . 

证 明 车 k > 1， 则 有 0 < Xo <k <L 由 于 当 0 < 入 < 和 时，f(A) > 0， 即 有 trJ > 
0. 又 因 adetJ > 0， 所 以 特征 方程 存在 正 实 部 的 特征 值 ， 故 此 时 (w*,v*) 不 稳定 ; 而 当 和 Xo < 
入 < ko Bt, fO) < 0， 即 有 trJ < 0， 此 时 特征 方程 存在 两 个 负 实 部 的 特征 值 ，(w*,v*) 局 部 
渐 近 稳定 . 

By = alà) ibA) 为 特征 方程 的 两 个 根 ， 则 有 





243 _ 2 = 
NE ec 9, (a) = 5/20 - at) - 


d2(\2k3 — 2Ak? + k — 1)? 
(1 — Ak}? 








经 计算 
a(Xo) =0,， aN, = -NR < 0. 


#242 Poincaré-Andronov-Hopf 分歧 定理 加 可 知 ， 当 和 = AoW, ABE (2) 在 平衡 点 (w*,v*) 处 产 
AE Hopf i. 
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AG=u-ut,T=v-v%. ATH, Alu, v KRU, T 那么 系统 (2) 可 变换 为 


SN 
UE v G? (u,v, A) 





其 中 
Gl(uwu A) = k(1 — Ap)au2 + Auv — k?(1 — Ak)tu? + O (luf, |u|* |u|), 
G?(u, v, A) = —uv + O (lult, ul |v). 
作 变 换 
uw 1 0 T 
v] \MN y |? 
ji A?k? — 2Ak? +k +1 1— Mk z 
M=- L, N=- VAD-T. 
A (4) 可 得 
Tt = a(A) —B(A) T n F(x,y, A) 
v: pa) aA J \y F*(z,y,d) 人 
其 中 


Fl(z, 和) = A202? + Ansy 十 42iz28 + 4soz3 十 Of(lz| 和 |z|?ly)， 


M 1 
F?(z,y, A) = -NF nÀ) a (Mz? + Nay) x 


Az = k(1 — àk)? +AM, Arr =AN, Agr =0, Ago = —k?(1 一 入 6) 
下 面 通过 判断 a(o) 的 符号 来 讨论 周期 解 的 方向 及 稳定 性 ， 由 文献 [10] 知 


1 
a(ào) = 75 [Face + Feyu + Fizy + Fi] 


1 1 
+ 16B00) e 
这 里 各 导数 均 为 导 函 数 在 (0, 0, Ao) 处 的 导数 值 . 根据 F (r y, rd), Fely A 的 表达 式 ， 显 然 有 
Fly = Fisy = Fyyy = Fyy = Foy = 0- 


注意 到 Mok = 1 一 k-3， 进 一 步 计 算得 


Fi, + Fi) = F2 (Fa, + F?) at ee ee Fyy Fy,l ’ 


Flaa = 6Agolaca, = —6k2(k-?)* = 一 6， 


T 


F} = Aiila=o a AoNo, Fe 一 一] 一 Ao Mo = 0, 
FL, = 2A20|y=n9 = 2(k $+ AoMo) = —2Aok, 


2M 2) yy 
F}, = Fy (A20 十 1)|>=xo = Ne ika, 
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其 中 
S 1— ok 
Mo = M|r=ro = 一 A0 3 No :一 NIx=xo Sa V Ao(1 = Aok). 





入 2 
HERA ao) 的 表达 式 ， 整 理 得 i 
ao) = 75 Flee + gogig FiFi- Fon] 
= -$ = Samy (oN = aoe?) 
z -$ z srk Ag A — dok)-#( - AJ? (L — Aok)? + 22 R(1 — Mok)3) 
2 -Š = 3 (QkX0 —1)Vk = n ZUR <0. 





注意 到 ww(Xo) < 0. 因此 由 Poincar6-Andronov-Hopf 分 歧 定理 知 ， 该 Hopf 分 歧 为 次 临界 方 
向 ， 且 周期 闭 轨 渐 近 稳 定 . 


3 ”扩散 系统 的 Hopf 分 歧 及 稳定 性 
在 常 微 分 系统 讨论 的 基础 上 ， 考 虑 带 扩 散 的 动力 系统 


Ut — divuzz = Mw — Ti, TE (0,7), t>0, 





Vt — dovVzz = 1 — uv, xré(0,7), t>0, (5) 


uz = Vz = 0, z=0,7, t>0. 


X = {(u,v) € H?(0,7) x H?(0, 7) : ue = vs =0, =0,7}, Xc=X +iX. 


系统 (5) 在 (ux, Va) 处 定义 在 Xc 上 的 线性 化 算 子 具有 形式 


TARA 


其 中 D= diag(d;, d2)» J 由 (3) 给 出 . 
我 们 知道 ， 特 征 值 问题 -oz = ud, x € (0,7), 各 = 0, z= 0 具有 一 序列 简单 特征 
ÍË Un =n? (n = 0,1,2,…)， 对 应 的 特征 函数 为 
1 
Vr’ 


2 
V2 osna), n>0. 


n=0, 


pn (7) = 


考虑 工 的 特征 方程 
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令 
$ 2 an 
= cos(nx), 
代入 特征 方程 得 
5 (Ja(à) — uI) ( ne ) cos(nz) = 0, 
n=0 n 
其 中 
= = à? 
Jn (A) := J — Wn D = | M0 TA 5 
inal -IF T dzin 
入 TAR 


A |Jn(A) — wl] =0, Wy LOA) 的 特征 方程 为 
-Tap + Dn =0, n=0,1,2,---, 


其 中 


MA2k3 — 24k? 十 大 一 1 
Th = Ae ew tk) — (dı + d2)ln, 


Dn = A(1 — àk) + didau? + rin — d2Ak(1 — Ak) ln. 

定理 2 设 0<k<1. #d, do Wd > Ak(1-— Ak) 3È di/dz2 > k(1—Ak)?, WABE (5) 的 
惟一 正常 数 解 (ux, Us) 渐 近 稳定 . 

证 明 只 需 说 明 对 任意 7 > 0, (A) < 0，Di(A) > 0. 当 0 < < 1 时 ,To() <0. A 
WET; (A) < To(A) <0. Æ di > Mk(1 一 和 Ak)， 则 对 任意 > 0， 有 


dA 
1— rR 





Di(A) = A(1 — Ak) + dzu; [diy — Ak(1 — AK)] + > 0. 
# dı/d2 >k(1—Ak)?, WHER > 0， 有 


入 


D;(A) = A(1 — Ak) + didau? + 
当 k > 1 时 ， 经 计算 知 
To(Ao) =0, Do(Xo) > 0， a(Ao)=0, oN0) <0, TM) < To(Ao) =0,j>0. 
注意 到 Xo = k7? (Vk —1). 车 di > k-1(VE 一 1) 或 di > do, W4I>OK, A 


D;(Ao) Ao(1 = àok) + dy dos + Nok? dip; = Nok? do; 


ll 


Ao(1 一 Aok) + dap; (dips; 一 Nok?) + Nok? dys; > 0, 


D;(ào) = Ao(1 — Aok) + didy? + Agk? dip; — Aok? dap; 
= No(1 — Aok) + didgp? + Aok? (di — de); > 0. 
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Att L(Ao) 除 一 对 共 轿 纯 虚 数 特征 值 外 ， 其 他 特征 值 均 具有 负 实 部 。 由 文献 [9] 中 的 Hopf 4 
歧 定 理 知 ， 系 统 (5) FE (Xo; us, Va) 处 产生 Hopf 分 歧 . 

定理 3 设 k>1. 若 di > k-1(Vk 一 1) 或 di > dz, MRA (5) Æ (ux, Vx, A0) 处 产生 Hopf 分 
歧 ， 且 该 Hopf 分 歧 为 次 临界 的 ， 周 期 闭 轨 渐 近 稳 定 ， 其 中 Xo = k-2 (Vk 一 1). 

证 明 利用 文献 [6,9] 中 的 方法 来 讨论 Hopf 分 歧 的 方向 与 稳定 性 . 

定义 L* (An) 是 线性 化 算 子 L(X0) 的 伴随 算 子 


e(z) 
v Vex v 


其 中 


令 


3 -3 PERE 
f% 1 | NY; veka | Ap KTE iN? 
1:5 F ara eRe: P 2 l ， 
bo —rAp +iAo 2k-4 bò i 


则 q, q* 满足 L(Ao)¢ = 180g, L*(Xo)q* aa —ifoq", (q*,@) = 1, (q*,@) = 0， 其 中 
t= B00) =K, (ast) = | bdr 
表示 L2(0,7) x L7(0,7) 中 的 内 积 . $ 


f(u,v, A) = àw 一 glu, v, A) =1—w. 


u 
1+ ku’ 
在 (Us, Ve, 0) 处 计算 得 

fuu = kU Ak) =2k-?, fao fos = 0, 

fuuu = —6k?(1 — Aok)* = —6, fuw=0, fuw = 0, fow = 9, 


Guu = 0, Qu = -1, Jw= 0, Quuu = Juuv = Juvv = gouu = 0. 


Qu = ‘ cos (nz), Qag = k |an Cagg = | Josten, 


Cn = fuua? 十 2 fuvdnbdn =F fuvb2, dn = guuQ2 二 29uvanbn + guo 了 2， 





定义 


其 中 


En = fuulanl? 十 fuv (anbn or Gnbn) + fuvlbnl?, 
fn = Gialeal? =F guv (anbn +F Gnbn) + guvlbn|?, 
gn = fuuulan| an F funy (2ļan| bn + až bn) F Juw (2|bn] an + b? Gn) H fovvlbn| bn, 


hn = Juuulan| an + Juuv (2ļan| bn + a2 bn) + guvv (2|bn|?an + ban) + Jvvv|bn| bn. 
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当 n = 二 0 时 ， 有 


2(k-? — 1+iN 2k-2) 2(k-4 — 1) -6 
Qu = E NE » Qa = 2 ， Caga = : 
2(A5 — ig 7&7 4) 29 0 


进一步 计算 得 (or, Cyog) = -3 + AZKE, B 


1 1 


ld, Qu) = k7? -1+ iN Sk 一 iĝ kā, (q", Qag) = k7} -1- ig ke, 
(T, Qua) = K7? -1 +iAJ Sk NR, (g, Qaa) = ho} —1 4 ing. 
定义 
pa : -1 x = 
q:= [2i6o7 = L(Ao)] w20, d :一 —[L(Ao)] Wii, 
其 中 
wo 一 Qaa Z (q*, Qad 一 (T, Quo)0， w11 = Qag Sx (q", Qaa)d z (J, Qaq)d- 


经 计算 ，wzo = 0, wi = 0， 从 而 有 (gq*, Qa) = (gq*, Qaa) = 0. 
由 文献 [9] 知 ，Hopf 分 歧 的 方向 由 Re(ci(Xo)) 的 符号 所 决定 ， 其 中 


i * * * 1 * 1 * 
cı (ào) = TAKI Qaa) ` (9%: Qaa) + (0°, Qaa) + 5 (9"s Qaa) + 5 (9" Caga). 
由 前 面 计算 得 
i 。 3 
Re 人 (ci(Xo)) 一 apap lia Qaa) - (q ,Qaa)] 一 了 
0 
= i(-inge i ri + sido PKT — sity k ag 3 
1 a ee Le r,t 
= SERS NT t+ 5A k 2 
= lkt lAlk-1(k3 — 
aa — hi + 5497 (KF - 1) 
Sneaky tet S 1) GS 
2 2 
1 lı 
一 二 了 <0 
由 文献 [6,9] 知 ， 车 


ggo (Xo)) < 0, 


则 Hopf 分 歧 方 向 是 超 临 界 的 ， 若 


1 
We) Re (er 0)) > 0, 
WW Hopf 分 歧 方 向 是 次 临界 的 ; HEE Lo) 其 他 所 有 特征 值 均 具有 人 负 实 部 ， 则 当 Re(c1(Ao)) < 
0 时 ， 周 期 解 渐 近 稳定 ， 当 Re(cl(Xo)) > 0 时 ， 周 期 解 不 稳定 . 注意 到 w(Xo) < 0， 因 此 
该 Hopf 分 歧 方向 是 次 临界 的 ， 且 由 该 Hopf 分 歧 产 生 的 周期 解 渐 近 稳定 . 
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注 1 对 于 扩散 系统 (5)， 我 们 只 是 得 到 了 空间 齐 次 的 Hopf 分 歧 周 期 解 (对 应 于 n = 0 的 模 
式 )， 该 系统 是 否 存在 非 齐 次 的 周期 解 有 待 于 进一步 的 研究 ， 事实 上 ， 空 间 齐 次 的 周期 解 与 常 
微分 系统 中 产生 的 周期 解 是 一 致 的 . 但 是 ， 对 于 扩散 系统 ，Hopf 分 歧 产 生 的 条 件 要 求 苛刻 ， 
周期 解 稳定 性 的 判断 也 更 为 困难 . 


4 ”数值 模拟 


这 一 节 ， 我 们 利用 Matlab 软件 给 出 数值 模拟 结果 ， 以 验证 定理 1 至 定理 3 中 的 结果 . 

常 微分 系统 包含 两 个 参数 : 入, k， 选 取 k = 4 > 1， 那 么 临界 值 Xo = 0.125， 由 定理 1 知 ， 
当 入 = Xo 时 ， 系 统 (2) 产 生 次 临界 的 Hopf 分 歧 ， 且 该 Hopf 分 歧 周 期 解 渐 近 稳定 ， 选 取 入 = 
0.12, WALA. 大 量 的 数值 模拟 表明 ， 当 入 > Aot, (us, ws) 整体 渐 近 稳定 ， 选 取 入 = 0.15, 
KALA. 

偏 微分 系统 包含 四 个 参数 : di, d2, A, k， 首 先 考虑 0 < k < 1 情形 . 选取 k = 0.8. HE 
H241, FF dı > 0.16 3 dı/d2 > 0.032， 则 (ux, v) YDET. 选取 入 = 1, dı = 0.2, dz = 10, 
见 图 2， 大量 的 数值 模拟 表明 ， 当 0 < k < 1 时 ， 系 统 (5) 不 会 产生 Hopf 分 歧 . FARM, 
当 qd1/ds 充分 小 时 ， 系 统 (2) 存在 稳定 的 非常 数 正 平衡 态 解 ， 选 取 入 = 1, di = 0.01, d2 = 10, 
见 图 3. 

其 次 ， 考 虑 k > 1 的 情形 ,选取 k = 4， 那 么 临界 值 Xo = 0.125. 由 定理 3 知 ， 若 心 > 
0.25 或 di > das WHA = 和 No 时 ， 系 统 (5) 产 生 次 临界 的 Hopf 分 歧 ， 且 该 Hopf 分 歧 周 期 解 浙 
近 稳 定 . 不 妨 选 取 d1 = 0.02, dg = 0.01. 若 取 入 = 0.12 < 和 Ao， 系统 (5) 存 在 稳定 的 空间 齐 
次 的 Hopf 分 歧 周 期 解 ， 见 图 4. ÆRA = 0.05 < 0.12 < Xo， 系 统 (5) 仍 存在 稳定 的 Hopf 分 
歧 周 期 解 ， 并 且 周 期 变 大 ， 见 图 5. RA = 0.23 > Xo， 大 量 的 数值 模拟 表明 (us, v) = 
(2.875, 0.348) 渐 近 稳定 ， 见 图 6. 并 且 我 们 发 现 ， 当 di /az 充分 小 ， 系 统 (5) 存在 空间 非 齐 次 的 
周期 解 ， 选 取 入 = 0.23 > ào, dl = 0.01, da = 10， 见 图 7. 
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图 1: 参数 k 二 4 > 1, 入 二 0.12， 系 统 (2) FER Hopf DAAMA (E); 
A= 0.15, (us, va) 渐 近 稳定 (BS), FHR ult) LYK v(t) 
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图 2: 参数 =0.8<1, \=1, dı =0.2>0.16, d2 = 10, 
SEH KP I (ux, Uw) = (5, 0.2) 渐 近 稳定 
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图 3: 参数 上 =0.8<1 \=1, dı =0.01, d2 = 10， 正 常数 平衡 解 (uvus ) 不 稳定 ， 
系统 (5) 存在 稳定 的 非常 数 正 平衡 解 





t-axis 0 0 x-axis 


图 4: 参数 k 二 4> 1, 入 =0.12 < 0.125, dı = 0.02, d2 = 0.01, 
系统 (5) 存在 稳定 的 Hopf 分 歧 周 期 解 
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图 5: 参数 k 二 4 > 1 入 = 0.05 < 0.12 < ào, di = 0.02, dz = 0.01, 
系统 (5) 存在 稳定 的 Hopf 分 歧 周 期 解 ， 且 相 比 图 4 周期 增 大 





图 6: 参数 k= 二 4>1, 入 = 0.23 > 0.125, di = 0.02, dz = 0.01, 
SEH KF ARIF (us, Ux) = (2.875, 0.348) 渐 近 稳定 





图 7: 参数 k= 二 4> 1, 入 = 0.23 > 0.125, dı = 0.01, dz = 10, 
系统 (5) 存在 不 稳定 的 空间 非 齐 次 的 Hopf 分 歧 周 期 解 
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Hopf Bifurcation and Stability for an Autocatalytic System 
MA Xiao-li 


(Department of Mathematics and Physics, Xi’an Technological University, Xi’an 710032) 


Abstract: This paper concerns a chemical model with an autocatalysis and saturation law. First, the 
subcritical Hopf bifurcation is obtained for the ordinary differential system and the induced periodic 
solutions are locally asymptotically stable. Then the diffusive model is considered. When the diffusive 
coefficients satisfy certain conditions, the subcritical Hopf bifurcation is also obtained and the spatially 
homogeneous periodic solutions are asymptotically stable. At last, numerical examples simulated with 
Matlab are shown to support and strengthen the analytical conclusions. 
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